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Abstract 

The perturbation expansion of the MS coupling in the SU(iV) gauge the- 
ory in powers of the bare lattice coupling is extended to two-loop order. A dis- 
cussion of our results has already been published elsewhere. We here describe 
the calculation in some detail. It relies on the background field technique and 
a recently introduced coordinate space method to evaluate lattice Feynman 
integrals. 
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1. Introduction 



Our aim in this paper is to compute the coefficient c?2( s ) m the expansion 
a Ms( s / a ) = a o + ^i(s) al + d 2 (s) al + . . . , (1.1) 

which relates the running coupling in the MS scheme of dimensional regular- 
ization to the bare lattice coupling cio in the pure SU(iV) gauge theory. The 
motivation for this calculation and references to earlier work were given in 
a previous publication [1]. We assume the reader is familiar with this paper 
and also refs.[2,3], where we discuss the background field formalism in lattice 
gauge theory and describe some coordinate space techniques to evaluate lattice 
Feynman diagrams. 

In sect. 2 the strategy of the computation is outlined. After summarizing 
some results for the dimensionally regularized theory in sect. 3, we write down 
the Feynman rules for the lattice theory in the presence of a background field 
(sect. 4). We then first calculate the one-loop diagrams contributing to the 
background field 2-point function (sect. 5). The computation of the two-loop 
diagrams is discussed in some detail in sect. 6. We finally collect our results 
and obtain c?2( s ) (sect. 7). 



2. Background field technique and renormalization constants 

The probably most effective way to compute the renormalization con- 
stants required to match different regularization schemes is based on the back- 
ground field technique. The method has previously been employed to calculate 
the one-loop coefficient di(s) [4-6] and its applicability to lattice gauge theory 
at all orders of perturbation theory has been established in ref.[2]. We have 
thus decided to make use of this technique to compute d,2(s). In the following 
the notation is as in ref.[2]. 

Let r( J ' fc -') be the bare vertex functions in the lattice theory and r^' fc '^ the 
renormalized vertex functions in the MS scheme of dimensional regularization 
at e = 0. The lattice vertex functions depend on the bare coupling go, the bare 
gauge parameter Ao and the lattice spacing a. In the case of the MS vertex 
functions the parameters are the renormalized coupling g, the renormalized 
gauge parameter A and the normalization mass ji. 
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As discussed in ref.[2] the vertex functions T^' k '^ and r^' fc '^ are related 

by 

T U,k,i) = Z k J 2 Zi T<- j > k >Q + 0(a), (2.1) 
where, for the parameters g and A on the left hand side, one should substitute 

g = Z^Z 3 3 /2 g , (2.2) 
A = Z 3 \ . (2.3) 

The renormalization constants Z\, Z% and Z% are formal power series in go 
with coefficients depending on Ao and a/i. 

The couplings appearing in eq.(l.l) are ay^(p) = g 2 /4ir and ao = g^/Ait. 
So if we define the coefficients Z^ through 



i^3- 3 = l + £2 ( V, (2-4) 
i=i 

it follows that 

d 1 {s) = -4vr Z {1) , (2.5) 

Ijl=s I a 

d 2 {s) = dtisf - (4vr) 2 Z^ . (2.6) 

/i=s/ a 

An important point to note here is that the coefficients Z^ (and hence the 
di(s)) are independent of the gauge parameter. 

To calculate Z^ and Z^ we consider the background field 2-point func- 
tion 

rk 2 '°' 0) (p, -p);1 = -^ a6 (<W 2 - - Mp)}/9 2 (2.7) 

and the inverse gluon propagator 

4°' 2 ' 0) (p, -p)f v = -S^US^p 2 - p, PlJ )[l - uj R (p)] + Xp, PlJ }. (2.8) 

The tensor structure of these functions, as given above, is implied by the 
symmetries of the theory. The Lorentz invariant amplitudes z^r(p) and ^r(p) 
are of order g 2 . On the lattice we define the analogous amplitudes v{p) and 
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oj(p) through 



E r(2 '°' 0) (p,-p) 



ab 
liji 



(2.9) 



E r(0 ' 2 ' 0) (p,-p) 



ab 
lili 



ry{3[l-^(p)] + A }, 



(2.10) 



where p^ = (2/ 'a) sin (ap^/2). They are invariant under lattice rotations and 
reflections, but not under continuous Lorentz transformations. 

In the rest of this section all contributions to z/(p) and oj(p) which vanish 
in the continuum limit a — ► are neglected. From eqs.(2.1)-(2.3) we then 
infer that 



Note that Z\ and Z 3 are only implicitly determined by these equations, since 
vr(p) and o>r(p) are functions of g and A, which in turn are given by eqs.(2.2) 
and (2.3). Introducing the loop expansions 



ZlZ^ = [l-v(p)]l[l-v K (p% 




Z 3 = [l-oj R (p)]l [l-uj(p)]. 




oo 





1=1 



oo 





1=1 



and similarly for the lattice amplitudes, 



oo 





1=1 



oo 





1=1 



it follows that 
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z^ = Lw(p)-JV(p)} _ (2-18) 

At two-loop order the renormalization of the gauge parameter must be taken 
into account and one obtains 

Z^ = U\ P )-^\ P ) + X d -^ziA . (2.19) 

I J A=A 

It is now clear that to compute c?2( s ) it suffices to work out o>r(p) and oj(p) to 
one-loop order and vr(p) and v(p) to two-loop order. While z/^ 1 (p) is required 
for general A, we may choose the Feynman gauge A = Ao = 1 in all other cases. 
The computations are significantly simplified by this choice, especially in the 
lattice theory. 



3. Computation of ojr and vr 



Starting from the total action, eq.(2.6) of ref.[2], it is straightforward 
to deduce the Feynman rules for the dimensionally regularized theory with 
background field. The relevant diagrams have been listed in refs.[7,8], for 
example. The computation of these diagrams is standard so that here we only 
quote the final results in the MS scheme [9,10]. 

At one-loop order the background field propagator is given by 

^ N f 11 , . 2 . 2 , 205 3 1 ] , , 

zA j (p) = <^ —l n (n 2 p 2 ) + + _ + — } , (3.1) 

R KFJ 16vr 2 \ 3 VP /F ' 36 2A 4A 2 J ' v ; 

while for the gluon self-energy one obtains 

oji 'ip) = <[ Mil 2 p 2 ) + — + — + — } . (3.2) 

R KFJ 16vr 2 \V 6 2A7 VP /F ' 36 2A 4A 2 J v ; 

vr{p) has previously been computed to two loops by Ellis [11]. His result 
was checked independently by van de Ven [12] and by one of the present 
authors. Prior to the calculation of Ellis, the divergent parts of the diagrams 
(from which the two-loop coefficient of the Callan-Symanzik /3-function can be 
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extracted) had been computed in refs.[7,8]. In the Feynman gauge the result 
is 

= (— 

A = l V 167T 2 



81n( M 2 /p 2 )-6C(3)+^}, 



f3.3) 



where ((z) denotes Riemann's zeta function (ref.[13], §9.522). 



4. Lattice perturbation theory 

In the remainder of this paper we discuss the computation of the lattice 
amplitudes v{p) and oj(p). We begin by describing the Feynman rules for the 
lattice theory, starting from the total action given in subsect. 5.1 of ref.[2]. 
For notational convenience we now set a = 1. 

4-1 Vertices 

The total action Stot of the lattice theory can be expanded in powers of its 
arguments according to 

S tot [B,q,c,c] = S Tip 

j,k=0 1=0,1 J ' 

x (27r) 4 ^( Pl + ... + S,) V^'Hpu- ■ • , s,)l\-;; d J k 

X B%(- Pl ) • • .B% (- Pj ) qll(- qi ) • ..qli (-?*) 

x l c H-n) . ..F'i-n) ~c dl (-si) ■ ■■c dl (-s l ) (4.1) 

(for unexplained notations see ref.[2]). The vertices V^' k ' 1 ^ are defined through 
this expansion and the requirement that they should be invariant under per- 
mutations of the momenta and indices associated with the same type of fields. 

A further set of vertices arises from the a priori measure in the functional 
integral. The jacobian appearing in eq.(5.10) of ref.[2] is equivalent to the 
"measure" action 

3 

Sm[q] = " E E tr i ln i J (9oq,(x))] } . (4.2) 
xeA 11=0 
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The associated vertices, Vm \ are obtained by expanding 

¥ M _f i r dv 

^M-L, kl J (27r) 4 • • • (27r) 4 



X (27r) 4 ^( Pl + ... + Pk ) Vi k \ Pl ,. . ., Pk )l\ 



a 1 ...a k 



Xq a a 1 A-Pl)---qZ(-Pk). (4.3) 



(k) ■ k 

Note that Vm is proportional to gfi and so is of higher order than the vertex 

y(o,fc,o)_ j n particular, the measure vertices do not contribute at tree-level of 

perturbation theory. 

For small values of j + k + 2l it is not difficult to obtain explicit expressions 

for the vertices V^' k ' 1 ^ by hand. For illustration the complete set of 2- and 3- 

point vertices is listed in appendix A. In our two-loop calculations vertices with 

up to 6 legs occur. The formulae for these are rather lengthy and it is advisable 

to generate them using an algebraic manipulation program. A possible way 

to proceed has been described some time ago in ref.[14]. Following these lines 

we have written MAPLE programs for all required vertices. The programs 

deliver an abstract algebraic expression in the colour tensors tr {T ai . . . T a ™}, 

the generalized Kronecker delta symbols 

{j _ J 1 if Mi = M2 = • • • = Mn, / 4 4 ^ 

\0 otherwise, 1 ' ' 

and the sines and cosines of the momenta. 
4-2 Graphical representation 

The Feynman rules for the vertex functions T^' k ' 1 ^ can now be set up in the 
usual way. It may be helpful at this point to consider the diagrams contributing 
to r( 2 ' ' ) at one- and two-loop order shown in fig. I in sect. 5 and fig. 2 in 
sect. 6. There are two types of internal lines, representing the propagation of 
the gluon field (wiggly lines) and the ghost fields c and c (directed dashed 
lines). The propagators associated with these lines are given by 



S^iS^ + iX- 1 -l)p,pjp 2 }/p 2 , (4.5) 
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— = S ab /p\ (4.6) 

P 

where the notation p^ = 2 sin^-p^) has been used (cf. appendix A). There are 
no factors associated with external gluon and ghost lines for which the obvious 
graphical notation is employed. 

Background field external lines are denoted by a wiggly line with a termi- 
nating cross. In particular, at tree-level the background field 2-point function 
is given by the diagram 



p,a,fj -p.b.v 



. v (2,o,o) {pj _ p)ab (4J) 



The vertices V^' k ' 1 ^ are represented by filled circles from which j + k + 21 lines 
emerge. The j background field lines are external and thus terminate with a 
cross. The notation is best illustrated by considering the 3-point vertices 

r,c,p 

= -^ (1,2,0) (p,?,rfc (4.8) 

p,o,p \ 

q,b,v 




V^ip^r)^. (4.9) 



The graphical symbol for the measure vertices Vm' is an open square as in 

p.oT'^.v = -^i 2) (P, ~pC- (4-10) 

In all these figures the momenta associated with the legs of the vertices are 
in-going momenta. Note that it is always minus the vertex yO'-M or vi k) 
which is associated with the corresponding graphical element. 

Without loss we may assign external and loop momenta in such a way 
that momentum is exactly conserved at each vertex (not just modulo 2ir). 
The loop momenta k\, . . . , k n must be integrated over the Brillouin zone with 
measure d 4 ki . . . d 4 k n /(2ir) 4n . Otherwise the rules are as in the continuum 
theory. 



/ r,c 



p.a.p \ 

v q.b 
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4-3 Computation of Feynman integrands 

At two-loop order the integrands associated with the diagrams contributing to 
the background field 2-point function may easily involve several hundred terms. 
To compute the integrands and to be able to manipulate them safely a fully 
automated approach is thus suggested. The algebraic calculation of Feynman 
integrands on a computer is standard in continuum theories. Similar programs 
may be used for lattice diagrams to work out the required contractions of 
colour and Lorentz indices. In particular, the contraction of the colour indices 
may be carried out by repeated application of the identities 

tr {T a XT a Y} = ^ tr {XY} - \ tr {X} tr {Y}, (4.11) 

tr {T a X} tr {T a Y} = ± tr {X} tr {Y} - \ tr {XY}, (4.12) 

which hold for arbitrary complex N X N matrices X and Y. 

Concerning the contraction of Lorentz indices we note that the lattice 
symmetries are not nearly as restrictive as the continuous Lorentz invariance. 
This is reflected by the presence of the generalized Kronecker delta symbols 
(4.4) in the algebraic expressions for the vertices. In general the result of 
the computation is hence not simply composed of scalar products of Lorentz 
vectors. One rather ends up with products of sines and cosines of various 
linear combinations of the momenta, with indices that are summed over or are 
equal to one of the external indices. To achieve a uniform representation of 
such terms it is sometimes useful to expand the products in sums of sines and 
cosines. The latter may then again be factorized into sums of products of 

k 2 , = 2 — 2 cos(A; |U ) and k^ = sin(A; |U ), (4-13) 

where k is any of the basic momenta. Since 

K = K- V k % ( 4 - 14 ) 

we may insist that there is at most one factor k^ in each product, for any 
given momentum k and index fj,. 

It should be quite clear from these remarks that the available programs to 
compute Feynman integrands in continuum theories are not directly applicable 
to the lattice theory. We have thus developed our own libraries of MAPLE 
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programs, which enable us to obtain the integrand associated with any given 
lattice Feynman diagram without manual calculations. One only needs to 
supply the diagram data in some appropriate form and to indicate which 
indices are to be contracted. Similar sets of programs have recently been 
described in refs.[15,16]. 



5. Computation of v^(p) and u^(p) 

At one-loop order there are 4 diagrams contributing to T^ 2 '°'°\p, —p) 1 ^ 
(see fig. 1). Accordingly we have 

i/^Cp) =J>(p), (5.1) 

Ao — 1 . 

8 = 1 

where the label i counts the diagrams. In the case of the "tadpole" diagrams, 
i = 3,4, the corresponding Feynman integrals are of the form 

while for i = 1, 2 the loop momentum may be chosen such that 

f I/ .( P ) = rj!iw k ± = k±h. (5.3) 

The functions Fi(p,k) are some complicated polynomials in the sines and 
cosines of the momenta. They are invariant under the lattice symmetries and 
we may assume that Fi(p, k) = Fi(p, —k) since the other factors in the integrals 
are invariant under reflections of k. 

We are interested in evaluating the diagrams in the continuum limit. Since 
the lattice spacing a has been set to 1, the limit amounts to taking p — ► 0. 
Our task then is to work out the asymptotic expansion of the integrals up 
to order p 2 . In the conventional language the terms of order p° represent 
a quadratically divergent contribution. These terms cancel in the sum of 
all diagrams. Terms of order p 1 (which would represent a linearly divergent 
contribution) are excluded by symmetry. In the continuum limit we are then 
left with the terms of order p 2 . 
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3 4 



Fig. 1. One-loop diagrams contributing to the background field 2- 
point function 

To compute the asymptotic form of the integrals we first expand the 
numerators Fi(p,k) of the integrands in powers of p. It suffices to keep all 
terms up to second order, since the remainder makes a contribution of order 
p 3 (or higher) to p 2 Vi(p) and so is negligible in the continuum limit. The 
numerators of the integrands then are of the form 

3 

f(k) + P^UAk) (5.4) 

with coefficients f(k) and /^(k) that transform like a scalar and a symmetric 
tensor under the lattice symmetries. 

In the case of the tadpole diagrams, i = 3,4, the symmetries allow us to 
substitute 

3 3 

E P^UAk)^\p 2 Y,f^(k)- (5-5) 

/i,i^=0 /i=0 

Inspection then shows that the integrals reduce to a linear combination of 1 
and the integral (B.l) introduced in appendix B. 

To evaluate the other diagrams, i = 1,2, we note that the integrals with 
integrands proportional to f(k) reduce to linear combinations of the integrals 
(B.10),(B.13) and (B.14). In the integrals with integrands proportional to 
fnis(k) we substitute 

f,Ak) = cS, 1J + flJk), (5.6) 

where c is a constant and /^(k) of order k 2 for k — ► 0. The terms proportional 
to c reduce to the integral A(p) [eq.(B.8)]. The remaining terms are sufficiently 
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regular that we can take the continuum limit of the integrals simply by setting 
p = in the denominator of the integrand. Taking the lattice symmetries into 
account these integrals reduce to a linear combination of 1 and the integrals 
(B.l) and (B.5). 

The steps outlined above are sufficiently systematic to be programmable. 
We have in fact written a set of programs which perform the necessary expan- 
sions and separation of terms. The integration, which amounts to matching 
the terms with the integrals listed in appendix B, has also been automated. 
The programs thus perform the complete computation. At no point does one 
need to copy and manipulate the terms by hand. 

Our final result then is 




where the constants Pi and P2 are defined in appendix B. 

The gluon self-energy ui^^p) has previously been calculated in ref.[17]. 
Apart from the diagram (4.10) which involves the measure vertex Vrfi , all 
diagrams contributing to the self-energy are obtained from those in fig. 1 
by replacing the background field external lines by gluon external lines. The 
computation of the corresponding Feynman integrals proceeds as in the case of 
the background field 2-point function. We only need the result in the Feynman 
gauge and there it is given by 

m/ n Af f 5 W 2 , 28) 1 

+ N {h p i + h + h} + ^ (5 - 81 

As an additional check of our programs we did in fact compute the full ampli- 
tudes r( 2 '°'°)(p, —p) a ^ v and r(°' 2 '°)(p, —p) a ^ v in an arbitrary gauge and verified 
that the calculated functions satisfied the Ward identities. 



11 



6. Computation of v^ 2 \p) 



We now describe the computation of the background field 2-point function 
at two-loop order in some detail. This is the most difficult part of our calcu- 
lations. We first outline the general strategy and then consider the diagrams 
one by one. The results are tabulated in subsect. 6.8. 

6.1 Preliminaries 

The diagrams contributing to r( 2 '°'°)(p, —p)"^ v at two-loop order are listed in 
fig. 2. In some cases the diagrams have been given a qualified label i.j, where 
i denotes the type of the diagram and j counts the diagrams of the given type. 
The total amplitude is then given by 

35 
i = 5 

where i>i{p) is the sum of the contributions of the diagrams of type i. 

In the continuum limit p — ► the amplitudes i>i{p) may be expanded in 
a series of the form 

3 P 4 f 

fviip) = c ,i + c lti E Hf + p ] C2 >« 

(6.2) 

The coefficients c n ^ have a simple dependence on N given by 

c n ,i = c%/N* + c™ + c™N*. (6.3) 

The second term in the expansion (6.2) is of order p 2 , but since it is not in- 
variant under continuous Lorentz transformations, we expect that these terms 
cancel in the sum of all diagrams. 

As already mentioned in sect. 4 the integrands of the diagrams are ob- 
tained by computer algebra. The general procedure to extract the coefficients 
c nt i from the corresponding Feynman integrals is as follows. We first expand 
the numerator of the integrand in powers of p. In most cases it is sufficient to 
work out the expansion up to second order, the remainder being negligible in 
the continuum limit. The terms are then classified according to their behaviour 
for small loop momenta. Terms proportional to sufficiently high powers of the 



z/( 2 )(p) 



ln(p 2 



+ C 3 i 



ln(p 2 



+ C 4 , 8 } + 0( P 4 ). 
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— > -- ' - — -<-.--' 

10.1 10.2 




Fig. 2. Two-loop diagrams contributing to the background field 2- 
point function 
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23.1 23.2 



24.1 24.2 



24.3 24.4 



Fig. 2. (continued) 
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loop momenta are easily dealt with. The expansion of the integral for p — ► 
may in these cases be obtained simply by expanding the integrand up to sec- 
ond order of p. One then ends up with a linear combination of integrals of the 
type 



where P(k, q, r) is a polynomial in the sines and cosines of the momenta k, q, r 
and n,m,l are positive integers. Such integrals can be computed straightfor- 
wardly using the coordinate space methods described in ref.[3]. 

The terms left after this must be treated individually and are eventually 
reduced to the special one- and two-loop integrals listed in appendix B and 
C and integrals of the type (6.4). As in the case of the one-loop amplitudes 
considered in sect. 5, all algebraic manipulations of the integrands are pro- 
grammable. The evaluation of the integrals of the type (6.4) (of which we had 
on the order of one hundred in our list) and those defined in appendix B and 
C is the only place where a numerical computation is required. 

In the following subsections a few further indications of how to proceed 
in the various cases are given. It is not our aim to describe all manipulations 
in full detail. We just discuss some of the key steps that have been taken, 
hoping that this will be helpful to anyone interested to compute two-loop 
lattice Feynman integrals. 

6.2 Factorizable diagrams (i = 25 — 29,34,35) 

These diagrams can be written as sums of products of one-loop integrals which 
can be computed following the lines of sect. 5. In the continuum limit all these 
integrals can be expressed analytically in terms of the constants Pi and Pi- 

6.3 Bigmac diagrams (i = 30 — 33) 

The integrals in this class are of the form 



where r = —k — q if i = 30, 31 and r = —k — q — p otherwise. The numerator 
Fi(p,k,q) is a polynomial in the sines and cosines of the momenta. 

For i = 30, 31 the expansion of the integrals in the continuum limit p — ► 
is obtained simply by expanding the numerator of the integrand up to sec- 
ond order in p. As a consequence of the lattice symmetries, the resulting 




r = — 



k-q, 



(6.4) 




(6.5) 
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expression, 

3 3 

f( k ,l)+^2p^f^ k ,l)+ ^2 P»Pvf»v(k,q), (6.6) 

/i=0 /i,i^=0 

may be replaced by 

3 

f(k,q)+\p 2 Y / f^(k,q)- (6.7) 

fj,=0 

One is then left with integrals of the type (6.4). 

Before the integration is performed it is worthwhile to simplify the inte- 
grand as far as possible. This may be achieved by rewriting the numerator of 
the integrand as a polynomial in k 2 , q 2 ,, k p and q a . Using the trigonometric 
identities (4.14) and 

Mm = W'l -K~ + IK€ ( 6 - 8 ) 

(with r = —k — q), all terms involving k p and q a can be eliminated. The 
resulting expression is a sum of monomials in k 2 , q 2 and f 2 . As it turns out 
there are just a few such terms left, and these are easily worked out using the 
coordinate space method. 

In the case of the diagrams no. i = 32,33, we expand Fi(p,k,q) up to 
order p 2 . Higher order terms are negligible since the scalar bigmac integral 
discussed in appendix C is absolutely convergent at p = 0. The numerator 
then is of the form (6.6). 

The terms of order p° may be integrated as follows. We first split the 
integral in three pieces by substituting 

f(k, q) = /(0, 0) + [f(k, 0) - /(0, 0)] + [f(k, q) - f(k, 0)]. (6.9) 

The first integral is proportional to the scalar bigmac integral (C.l). The 
second term is of order k 2 for small k. In this part of the integral we substitute 
k — ► k — p (so that r — ► —k — q) and then obtain the desired expansion of 
the integral for p — ► by expanding the integrand to second order. The third 
part of the integral is treated in the same way, the substitution here being 
q — ► q — p. In all cases one ends up with integrals of the type (6.4). As above 
the lattice symmetries and trigonometric identities should be used to simplify 
the integrands. 

The integrals associated with the order p 1 and p 2 contributions to the 
numerator (6.6) are easy to work out. We simply insert the expansion of 1/f 2 
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to order p 1 and p° , respectively. Again one is then left with integrals of the 
type (6.4). 

6.4 Tadpole diagrams (i = 12 — 18) 

These diagrams consist of a gluon or ghost loop with a self-energy insertion. 
The corresponding contributions i>i{p) are formally given by 

fMp) = j\wf {P) ~ 2HP ^ (6 - 10) 

where k is the momentum flowing through the loop. For i = 16 the function 
Fi(p,k) is a polynomial in p 2 ^ and k 2 . In all other cases it is an integral over 
a second loop momentum q. 

From the structure of the ghost-gluon vertices it follows that Fn(p,k) 
and Fis(p,k) are of order k 2 at k = 0. The integral (6.10) is thus well-defined 
in these cases. The same is not true for the other diagrams, where instead one 
may show that 

F t (p,k) = r t F w (p,k) + 0(k 2 ) (6.11) 

for some constants r^. Since the sum of all self-energy insertions is of order k 2 
we have 

16 

r t = 0. (6.12) 

i = 12 

The conclusion then is that while the integrals (6.10) do not separately exist 
for i = 12 — 16, their sum is perfectly well-defined. Taking eq.(6.12) into 
account, it is clear that the same result is obtained if we define i>i{p) through 
the subtracted integral 

P 2 <P)= I S^C^r'i^^-r^eip.k)]. (6.13) 

The expansion coefficients c n ^ tabulated in subsect. 6.8 refer to this definition. 

Since the denominators of the integrands do not depend on p, the ex- 
pansion of the integrals in the continuum limit is obtained by expanding the 
numerators up to order p 2 . The resulting integrals are of the type (6.4), or a 
subtracted variant of these integrals, such as 




(6.14) 
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The coordinate space technique applies to integrals of this type as well. 



6.5 Eye diagrams (i = 22 — 24) 

The contributions from these diagrams are given by 



~2 r \ _ I d 4 £; d 4 g Fi(p,k,q) 



^ ) = /_^(2^I^^ (6 - 15) 

where s = k + p and r = —k — q. The numerator Fi(p, k, q) is a polynomial in 
the sines and cosines of the momenta. As usual we first expand Fi(p,k,q) up 
to second order in p, rewrite the resulting expression as a polynomial in p^, 
k 2 etc., and eliminate products involving k^ and q^ as far as possible. 

Using trigonometric identities and reexpanding in powers of p, the terms 
of lowest order in the momenta can be written as a linear combination of p 2 , 
k 2 , s 2 , q 2 and f 2 plus terms of higher order. The integrals associated with the 
second order terms in the linear combination reduce to the scalar eye integral 
discussed in appendix C or to integrals already treated before. 

The remaining terms are of fourth (or higher) order in the momenta. We 
first consider the terms of order p° . They come in the form of a polynomial 
f(k,q,r) in A; 2 , q 2 and f 2 . The polynomial is invariant under the lattice 
symmetries and is of order k 2 for small k. If we introduce the function 

1 r d 4 g f(k,q,r) 

the integral to be computed may be written as 

C d 4 k g(k) 



g{k) is smooth for all k that do not vanish modulo 2ir. Moreover it follows 
from the properties of f(k, q, r) that 

g(k) = g(0) + O(k 2 ). (6.18) 

It is now straightforward to show that 

g(k) t 2 r d 4 k ^ ktg(k) 



kv) = a + v) f — ^ - v r -—- V 
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"^(0) + O(p 3 ), (6.19) 

where use has been made of the lattice symmetries to simplify the expressions 
on the right hand side. The constant g(0) is a one-loop integral while the 
other terms appearing in eq.(6.19) are integrals of the type (6.4). 

The terms of first and second order in p are treated similarly. Some of 
these terms may not involve any factor of k 2 ^. In these cases the function g{k) 
is defined with a subtraction at k = 0. One then ends up with subtracted 
integrals of the type (6.4), while the contribution from the subtraction term 
reduces to a product of one-loop integrals. 

6.6 Ring diagrams (i = 5 — 11) 

The integrals associated with these diagrams are of the form 

f v d 4 k 1 

p 2 Vl (p)= / - — _ _ Fi(p,k), s = k + p. (6.20) 

AL > J_„ (2vr) 4 (fc2)2j2 h ^ v > 

For i = 9 the function Fi(p,k) is a polynomial in the sines and cosines of 
the momenta, while in all other cases it is an integral over a second loop 
momentum q. A subtraction must be applied to some of the integrals (6.20) 
to make them converge at k = (cf. subsect. 6.4). Explicitly, for i = 5 — 9, 
eq.(6.20) is replaced by 

P 2 Mp)= P ^-^—[F^k)-r % F,(p,k)l (6.21) 
J- v (Zirp (k 2 ) 2 s 2 

where the constants are determined through 

F t (p,k) = r t F 9 (p,k) + 0(k 2 ). (6.22) 

The sum of the subtracted integrals is the same as the sum of the unsubtracted 
integrals. 

The integrals no. i = 7 — 9, 11 factorize into products of one-loop integrals 
and so are easy to compute. The other diagrams may be worked out essentially 
following the steps taken in the case of the eye diagrams (subsect. 6.5). Some 
of the more difficult integrands are those with numerators of order p 2 having 
the lowest possible degree in the loop momenta. These terms can be reduced 
to the integrals E^ip) and S flll (p) discussed in appendix C. 
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6.7 Diamond diagrams (i = 19 — 21) 

The loop momenta in these diagrams may be assigned such that the corre- 
sponding Feynman integrals are given by 

2 , , r d4fc A Fi(p,k,q) 
fvdp)= / - — -t—-{t - } ' , (6.23) 

where the numerator Fi(p, k,q)is a polynomial in the sines and cosines of the 
momenta and 

k± = k ± g± = q ± r = —k — q. (6.24) 

Considering the symmetries of the denominator of the integrand, we may 
assume that Fi(p, k, q) is invariant under the transformation k,q — ► —k, —q. 

In the continuum limit p — ► 0, the scalar diamond integral is power- 
counting convergent and one may show that [18] 

* d 4 fc d 4 g 1 6 2 

(2vr) 4 (2vr) 4 k\PPq\q 2 - (^) 4 

A simple consequence of this observation is that the remainder in the expansion 

3 3 

Fi(p,k,q) = f(k,q)+ ^ P»Pvfv,u{k,q)+ ^ P t ,PvP P P<jf^ P <j(k,q) + 0(p 6 ) 

fi,iy=0 fi,iy,p, <7 = 

(6.25) 

makes a negligible contribution to the integrals (6.23) and so may be dropped. 
Moreover it follows that the contribution of the fourth order term is given by 

6 3 

7^4 C(3) J] P^PpP-W(0,0)/p 2 + O(p 4 ). (6.26) 

We are thus left with the first two terms in eq.(6.25). 

The terms of order p 2 are first classified according to their behaviour for 
small k,q. Inspection shows that all terms are of second (or higher) order. 
Those which are of second degree may be reduced to integrals considered 
previously by applying simple trigonometric identities such as eq.(B.15) and 

pk = \{k\ - k 2 _). (6.27) 
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The terms of fourth and higher degree can be expanded straightforwardly and 
reduce to integrals of the type (6.4) or subtracted integrals of this type. 
It turns out that the order p° terms may all be written in the form 

f(k,q) = (kq)g(k,q,r), (6.28) 

where g(k,q,r) is a polynomial in k 2 , q 2 and f 2 . The corresponding integral 
reduces to a product of one-loop integrals if g(k,q,r) is proportional to f 2 . 
When g(k,q,r) is a polynomial of lower degree 3 (or more) in both k and q, 
the expansion of the integral is obtained by expanding the denominator of the 
integrand up to second order in p. 

Up to perhaps a substitution k,q — ► q,k (which is a symmetry of the 
denominator of the integrand), the terms g(k, q, r) left over then are 

3 

k\ (k 2 ) 2 and Y,K f i- ( 6 - 29 ) 

fj,=0 

In the first two cases the integral may be reduced to integrals discussed before 
by inserting the identity (B.15). In the last case we expand the denominator 
after subtracting the integral over k at p = q = and end up with (subtracted) 
integrals of the type (6.4). 

6.8 Results 

We now list the coefficients c n ^ appearing in the asymptotic expansion (6.2). 
The coefficients c^^ and c^^ of the logarithms can be given analytically in 
terms of the constants Pi and Pi introduced in appendix B (table 1). The 
coefficients vanish for all i except for i = 5, 7, where we obtain 

N 2 f 1\ , 

As anticipated in subsect. 6.1 the associated terms in the expansion (6.2) 
cancel after summing over all diagrams. 

It remains to tabulate the coefficients co,j and c^. We begin by noting 
that the only non-zero contributions proportional to 1/N 2 to these coefficients 
are 

c (o) _ (o) _ _ 2 (0) _ j_ (6 sn 

c 4,14 — c 4,32 — ^ c 4,34 ~ 64 - ^U.Oij 
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Table 1. Coefficients C2,i and cs t i 



i 








5 


15 
4 


85 p 
36 P 1 " 


15 p 10 5 

2 2 9tt 2 32 


6 


5 

12 


Y2 P i ~ 


5 p 25 

6 2 144tt 2 


7 





35 p 


32 T 8 iv 


10 


1 

6 


2T P 1 " 


1 p 1 

3^ 2 12tT 2 


11 





_ T8 P 1 




19 


27 

8 


125 p 
144 ^1 " 


27 p 1213 
4 -* 2 11527T 2 


20 


1 

4 


23 p _i 
288 1 1 


1 p 1 43 

2 ^ 2 576vr 2 


21 


1 

24 


^6 P 1 " 


1 p 5 

12 2 128tt 2 


22 


9 
8 


119 p _, 


9 p , 9 , 7 
4 J 2 "T 128tt 2 ~T 32 


23 


5 

_ 24 




7 

1287T 2 


24 


1 

6 


-2T P i + 


1 p , 3 

3^ 2 + 32tt 2 


25 


-6 


"f P l + 


12P 2 + 2^2 +W4 


27 





91 p i 
24 P l + 


JL _ ijy-2 
64 8 


29 





T8 P 1 




32 





3 

16tt 2 




33 





1 

48vr 2 





c 2,i = c 3,i = for all i not appearing in the first column 



[cf. eq.(6.3)]. The contributions proportional to N° are listed in table 2. The 
constants Q\ and Q2 appearing there have been introduced in sect. 3 of ref.[3]. 
Their numerical values are 

Qi = 0.0423063684(1), (6.32) 
Q 2 = 0.054623978180(1). (6.33) 

In table 3 we finally list the contributions proportional to N 2 . Here only nu- 
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Table 2. Coefficients c ■ and c\ ■ 



7 


-Pi 


+ 1 


5 p 

96^1 ~ 


5 D 11 

8^2 1927r 2 




12 







1 p2 
32^1 


" t§2- p i + + hQt 


1 

512 


13 







_ A p2 i 
96^1 ^ 


~ 96 Pl ~ 96^! 




14 




l 

16 


1 p2 
24 ^1 


1 p 3 

8 ri 256 




27 


Pi - 


1 

8 


192 1 


f 8 P 2 + 96jr 2 




30 







3 p2 




1 

768 


32 







1 

192 






34 








13 p , 1 

96 1 "r" 512 





c o ! 4 ! ^ or 1 n °t a PP ear i n g in the first column 



merical values are quoted, because the analytic expressions for the coefficients 
in terms of a basis of integrals of the type (6.4) would be too lengthy. 

To guarantee the correctness of our results, all diagrams have been eval- 
uated by both authors, using independent sets of programs. This includes 
the computation of the Feynman integrands, the expansion of the associated 
integrals for p — ► and, finally, the numerical calculation of the integrals of 
the type (6.4). 

As a further check we note that the sum of the quadratically divergent 
terms vanishes, 

35 



J> Oj8 = 0, (6.34) 



as it should be. This can be verified from tables 2 and 3, but we have actu- 
ally been able to establish the cancellation of these terms algebraically. The 
coefficients of the logarithms also add up to the expected values. Our result 
for the amplitude v^ 2 \p) then is 



N 2 

ln(p 2 ) + co/iV 2 + ci + c 2 N 2 + 0(p 2 ), (6.35) 



A = l 327T 
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Table 3. Coefficients c~^] and c*^\ 



i 


(2) / , \9 

4;i x ( 47r ) 2 


eg X (47T) 2 


5 


-0.04037441(4) 


0.02348922(2 


6 


0.12388806(3) 


0.07287386(1 


7 


4.59937175(1) 


1.97571361(1 


10 


-0.09574219(1) 


0.01916905(1 


11 


-0.32549037(1) 


0.02049418(1 


12 


-0.09963952(2) 


-0.07789705(1 


13 


-0.07850005(1) 


0.03280987(1 


14 


-2.29968587(1) 


2.30107035(1 


17 


0.08108773(1) 





18 


0.16274518(1) 





1 Q 


94305394f 1 


1 7668960(9 


2U 


— 0.00054087(1) 


n nnooonA 1 ) / 1 

— 0.00823993(1 


21 


-0.02738405(1) 


0.00751124(1 


22 


-0.20607863(2) 


0.42631335(4 


23 


-0.09023515(1) 


-0.05439543(1 


24 


08461788f 1) 


-0 02030710f 1 


25 


0.07074594(1) 


-0.37533511(1 


27 


-4.74086363(1) 


-1.91757176(1 


29 


0.32549037(1) 


-0.02049418(1 


30 


-2.16944115(1) 


-0.09039338(1 


31 


0.23298215(1) 





32 


2.27248046(1) 


0.11846333(1 


33 


-0.23017351(1) 


0.00658707(1 


34 


2.37043182(1) 


-1.44106163(1 


35 


-0.16274518(1) 






c o ! 4 ! ^ or 1 no ^ a PP ear i n g in the first column 
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where the constants cj = X}?=5 c 4*l are gi yen by 

c = 3/128, (6.36) 
ci = -0.016544619540(4), (6.37) 
c 2 = 0.0074438722(2). (6.38) 

Note that to obtain the numerical values above, the algebraic results for the 
diagrams were summed before numerical evaluation. 



7. Final results 



where 



Combining eqs.(2.5),(2.6) and (2.17)-(2.19) with our results for the back- 
ground field 2-point function and the gluon propagator, we finally obtain 

UN 7T 

Ms) = - 1 —Hs)-^ + k 1 N, (7.1) 
d 2 (s) = d^sf - -j^H-s) + — + k 2 + k 3 N 2 , (7.2) 

h = 2.135730074078457(2), (7.3) 
k 2 = -2.8626215972(6), (7.4) 
k 3 = 1.24911585(3). (7.5) 

The coefficients of ln(s) in eqs.(7.1) and (7.2) are of course directly related 
to the first two universal coefficients of the Callan-Symanzik /^-function. The 
one-loop result was first obtained by A. and P. Hasenfratz [17]. It was subse- 
quently reproduced using the background field technique in refs.[4-6]. 

Finally we note that the constant k 2 appearing above also enters the 
expansion 

p = 1 - i r™» + ( jv2 - 1| {^ + ir + 8^) a ° + --- (7 - 61 
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of the (normalized) plaquette expectation value P. 

We are indebted to Anton van de Ven for performing various checks on our 
calculations, particularly of the results quoted in sect. 3 for the dimensionally 
regularized theory. We also thank Mrs. R. Heininger for producing the figures. 



Appendix A 

We here list the non- vanishing 2- and 3-point vertices defined in sub- 
sect. 4.1. The function P^ip) appearing in the expressions below is given 

by 

where p^ = 2sin(p |U /2). 

We first write down all 2-point vertices. The inverse of the gluon and 
ghost field propagators are 

^ ( °' 2 ' 0) (P, -P)t = ^ i P M + X °PM > ( A - 2 ) 

V^ '°' 1) (p,-p) ab = S ab p 2 . (A.3) 
For the vertices with background field external legs one finds 

V^(p,-pC = g^S ab P^(p), (A.4) 

V^ 1 ' 0) (p,- P )li = go 1 S ab PM, (A.5) 

and the 2-point "measure" vertex is given by 

Vi 2) (p,-p)£ = ±g 2 N6 a % v . (A.6) 

This completes the list of 2-point vertices. 

We now tabulate the 3-point vertices with in-going momenta p, q, r such 
that p + q + r = 0. For the 3-point gluon vertex one obtains 

V^' 3 '°Hp,q,r)l% = tgof abc {S,Ui^q) P cos(r,/2) 

+ f>v P (q - r),, cos(pj,/2) + Sp^r - p) v cos(g p /2)}, (A. 7) 
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while the ghost-gluon vertex assumes the form 

T/(°' 1 ' 1 )(p, q, r)l bc = zg f abc q, cos(rv/2). 



(A.8) 



In these equations f abc denotes the SU(iV) structure constants as defined in 
appendix A of ref.[2]. The remaining 3-point vertices are 



^ (3 ' ' 0) (p,?,r)^ = ^ 3 ^ ( ° l3l0) (P,?,r)- 



v^ Hp, q ,rC= g ^v^H P , q ,r); 



livp 



V^Hp, q ,r)l 



livp 



(A.9) 



+ Wr c {^ p pm - k p pm} , (a.io) 

+ \f abc {S, P Pv P (q)-S^P» P (r)} 

+ iX f abc {S, p q„e-^/ 2 - S^fpe-^ 2 } , (A.ll) 



(A.12) 



Appendix B 



The integrals listed below can be expressed in terms of two constants, 



P 2 = m - — -ln(m 2 )+ / - — — —„ }. (B.2) 

m ^0|(47T) 2 1 1 (27T) 4 (fc2+ m 2)2 J 1 J 

Their numerical values 

Pi = 0.15493339023106021(1), (B.3) 

P 2 = 0.02401318111946489(1), (B.4) 
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have been determined in ref.[3]. 

In the course of our computations a small number of one-loop integrals 
with zero external momentum occurred. Besides the basic integrals (B.l) and 
(B.2) the list of these integrals includes 




To establish eqs.(B.5)-(B.7) one makes use of the identity (4.14) and partial 
integration. 

In the case of the integrals depending on the external momentum p we 
are only interested in their asymptotic form in the continuum limit p — ► 0. 
The scalar self-energy diagram 

f v d 4 k 1 

A(p)= / - — -r^—, k± = k±±p, (B.8) 
U ) J_ v (2vr) 4 k\k 2 _ ±2^ \ ) 

has been worked out in ref.[3] using a momentum space subtraction method, 
the result being 

A(p) = j^[-Hp 2 ) + 2] +P 2 + 0(p 2 ). (B.9) 

Here and below the error term 0(p n ) stands for a remainder R(p) such that 
linip^o R(p)/\p\ n ~ e = for all e > 0. 

The momentum space subtraction technique may also be applied to the 
more complicated integral 

c ' Ai,) --[^ Ki \w_-w] iBM) 
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[cf. eq.(4.13)]. The calculations can be simplified by noting that 

3 

Y / p,c,Ap) = o, (B.ll) 

fj,=0 

and the result then is 

CM = (<W 2 " P»P») \-^Mp) - T^r-r 2 + ^Pi ) + 0(p 4 ). (B.12) 



12 yrj 18(47r) 2 72' 



To deduce the expansions 



v c] 4 k k 2 ( 1 1 

7^4t = Pi + p 2 \ --A(p) + —Pi <> + 0(p 4 ), (B.13) 

(2vr) 4 £^2_ ^ \ 4 W ^ 32 ' ^ U h y J 

v d 4 k (k 2 ) 2 



l + 0(p 4 ), (B.14) 



we insert the identity 



k 2 = 1(^2 + p ) _ i p 2 + ^pIK + 0(p 4 ) (B.15) 

fj,=0 

in the numerator of the integrands and expand in powers of p after cancelling 
factors of k 2 ±. 

We finally quote our result for three further one-loop integrals which were 
needed to evaluate some two-loop diagrams. They are given by 

v r] 4 k k ( 1 ^ ^ 

" -pA-^A(p) + ±:P 1 +0(p 3 ), (B.16) 



w (2vr) 4 k 2 s 2 ^ I 2 X1J 16 

v d 4 k k^k v _ 1 

7r (2vr) 4 (jfc2)2j 2 " 2(4vr) 2 p 2 

+ V {^-^j+O^ 2 ), (B.18) 



where s = k + p. 
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Appendix C 



The two-loop integrals listed here have been worked out using a mixture 
of momentum and coordinate space techniques of the sort discussed in sect. 4 
of ref.[3]. There are four basic integrals that must be evaluated in this way 
and thus require special attention. The simplest of them, already considered 
in ref.[3], is the scalar "bigmac" integral 



Dip) 



v d 4 k d 4 q 1 



_ v (2vr) 4 (2tt) 4 £2^2' 
In the continuum limit p — ► the asymptotic expansion 



k — q — p. 



(C.l) 



D(p) = A +p 2 



2(4ttY 



ln(p 2 ) + A 1 



0(p 4 



(C.2) 



holds, where 



A = 0.0040430548122(3), 
A x = -0.00007447695(1). 

The integral corresponding to the scalar "eye" diagram, 

^ d 4 k 1 



(C.3) 
(C.4) 



E(p) 



* (2vr) 4 jfc2j 2 



A(k), 



k + p, 



(C.5) 



has also been worked out in ref.[3] (the function A(k) is discussed in appendix 
B). In the continuum limit it is given by 



E(p) 



1 



Hp 2 )}' 



1 



2(4vr) 4 L v ^ /J (4vr) 2 
where, for the constant X\, the value 



47T 



Mp 2 ) + X 1 + 0(p 2 ), (C.6) 



Xt = 0.001123187953(1) 



(C.7) 



has been obtained. 
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Two further integrals, related to the "ring" diagrams, are defined through 



EM 

£M) 



d_ k k ^ k u 



n (27r) 4 (jfc2)2j2 
W d 4 k 1 



7r (2vr) 4 (jfc2)2j2 



s = A;+p, (C.8) 
^(l^cos^^C^) (C.9) 



[cf. eq.(4.13)]. Note that the function C ^{k) (which was introduced in ap- 
pendix B) is not invariant under shifts of the components of k by integer 
multiples of 2ir, but the integrand in eq.(C.9) is. The asymptotic expressions 
for these integrals have the tensor structure 



EM = Ei(p) + ^ E 2 (p) + 0(p 2 ), 

SM = hv M + ^ M + o(p 2 ). 



(CIO) 
(C.ll) 



For the scalar amplitudes E\(p) and £\{p) one finds 



Ei(p) 



M 



i i 



4 2(4vr) 4 



Hp 2 )]' 



(4tt) 



P 2 + 



2(4vr) 



ln(p 2 ) + X 2 



(47T) 



6 1 6(4vr) 2 



ln(p 2 ) + X 3 



with the constants X 2 and X3 being given by 

X 2 = 0.000583430648(4), 
X 3 = 0.000403986926(4). 



(C.12) 



(C.13) 



(C.14) 
(C.15) 
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For the other amplitudes, ^(p) an( i ^2(^)7 one obtains 



E 2 (p) 

Up) 



i 



2(4vr) 4 
1 

24(4vr)' 



ln(p 2 )- 

-MP 2 ; 



2(4vr) 



13 



24(4vr) 



P 2 - -P x 
6 



(C.16) 
(C.17) 
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